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Abstract 

We consider current observational constraints on the electromagnetic charge of dark matter. The 
velocity dependence of the scattering cross-section through the photon gives rise to qualitatively 
different constraints than standard dark matter scattering through massive force carriers. In 
particular, recombination epoch observations of dark matter density perturbations require that 
e, the ratio of the dark matter to electronic charge, is less than 10 -6 for mx = 1 GeV, rising 
to e < 10~ 4 for mx = 10 TeV. Though naively one would expect that dark matter carrying a 
charge well below this constraint could still give rise to large scattering in current direct detection 
experiments, we show that charged dark matter particles that could be detected with upcoming 
experiments are expected to be evacuated from the Galactic disk by the Galactic magnetic fields 
and supernova shock waves, and hence will not give rise to a signal. Thus dark matter with a 
small charge is likely not a source of a signal in current or upcoming dark matter direct detection 
experiments. 
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I. INTRODUCTION 



The nature of the dark matter (DM) remains a mystery. For DM in the MeV to TeV 
range, a wide variety of probes constrain the DM to be a Weakly Interacting Massive Particle 
(WIMP) which interacts with ordinary matter through suppressed couplings. These probes 
include direct detection of DM through nuclear recoils in underground detectors as well as 
indirect detection through DM annihilation to SM states in the sun (to neutrinos), in the 
Galactic center (to photons), and in the Galactic neighborhood (to charged particles). There 
are also significant constraints on DM couplings to ordinary matter through production and 
escape as missing energy at colliders. For a review, see 

Many of the most popular DM candidates naturally meet these stringent requirements. 
The neutralino from supersymmetry, for example, carries no electric charge and can interact 
only sub-weakly, via the Higgs or through small couplings to the Z boson, evading the most 
stringent constraints from LEP, Tevatron, and direct detection experiments such as CDMS 
j^J and XENON10 jsj. Its thermal annihilation cross-section is below the bounds for indirect 
detection through neutrinos, photons, or charged cosmic rays. While WIMP DM has escaped 
direct and indirect detection thus far, it may be within reach. Direct detection experiments 
are scaling up, the reach of the LHC will begin to encompass weak scale DM candidates 
soon, and Fermi will continue to constrain DM annihilation in dwarf galaxies, the Galactic 
center, and in the halo. 

At the same time, it is desirable to take as model-independent an approach as possible 
when constraining the nature of the DM. While in most popular models the DM carries no 
electromagnetic charge, periodically the notion of CHArged Massive Particle (a CHAMP) 
has reappeared in the literature jj-lSj- In some of the earliest discussions of CHAMPs, the 
DM carried a full unit of charge, but it was realized that this runs into a wide range of very 
stringent constraints from searches for heavy hydrogen to direct detection in underground 
labs. Some of these constraints may not apply if the CHAMPs are expelled from the disk via 
shock waves from supernova remnants and screened from re-entry by the Galactic magnetic 
fields jsj]. More recently, the possibility that DM carries a fractional or epsilon-charge has 



been considered and constrained via the CMB acoustic peaks 10J. Radio observations also 
constrain the electronic charge of the dark matter 111. In addition, the notion that the DM 



carries a "dark charge" has recently been considered 12Hl7l|. In these latter models the DM 
does not couple to the photon, but to a massless gauge force in the hidden sector. 

In light of the current understanding of structure formation and cosmological history, 
we determine how large the DM charge can be while remaining consistent with current 
constraints. We also consider direct detection signals from epsilon- charged DM, and 



determine whether it is possible to give rise to the signals in DAMA 18] and CoGeNT 



191 ] as discussed recently in 20|. Because we are answering a general question about the 



coupling of DM to the photon, we leave our discussion of models to a minimum. We note that 
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the discussion encompassed by this paper does bring to light a number of constraints that 
strongly disfavor some recent models in the literature. We comment on these models below 
where relevant. DM may also have a magnetic or electric dipole; this has been thoroughly 
considered recently [2l|, and we do not discuss it here. 

The outline of this paper is as follows. We begin with a brief discussion of models and the 
implications of this study for the viability of these models. We then review the relic density 
calculation before turning to constraints. We discuss halo shape constraints and the bound 
from scattering at recombination times. We discuss direct detection of charged particles in 
light of the signals from CoGeNT and DAMA, and the implications of the bounds discussed 
here for these experiments and models designed to fit them. Finally, we conclude. 



II. MODELS AND GENERAL CONSIDERATIONS 



Since DM that carries an electric charge must conserve £7(1)em, it must be a Dirac 
particle. There are a number of models in the literature where the DM carries a fractional or 
epsilon-charge. If a dark photon is massive and kinetically mixes with the photon, an epsilon- 



charge arises in Stueckelberg models [22[ on account of the unique form of Stueckelberg mass 
term. If, on the other hand, the dark photon is massless, kinetic mixing between the dark 
and visible photons induces an electric charge for the DM (or equivalently, a dark charge 
for visible states) [23[ . This mechanism is utilized for example in the Mirror Charged DM 
model proposed by [20| to generate the signals in CoGeNT and DAMA. We will see that 
the constraints we discuss here strongly disfavor such a model as the explanation for these 
signals. In either case, we denote the charge of the DM as ee. 

When determining the constraints on the DM charge, the essential features will be the 
irreducible coupling to the photon (and charged SM particles), and, more importantly, the 
velocity dependence of the scattering cross-section. For example, the Rutherford Scattering 
cross-section of DM off DM through a photon is 

da XX «em e4 /,N 



da " m^ r 4 el sin 4 (0*/2)' v ; 

where m>x is the DM mass, v re i is the DM relative velocity, and 9* is the scattering angle in 
the center-of-mass frame. Likewise, the scattering cross-section of DM off baryon is 



da X b «g m e 2 



da 4/^ r 4 el sin 4 (^/2)- 
where //& is the DM-baryon reduced mass. 



(2) 



The important point phenomenologically is the very large enhancement in the scattering 
cross-section at low velocity, giving a hint for where to look for strong constraints on 
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DM charge. Galactic constraints, where the DM has been heated through collapse and 
virialization, as we will see, tend to give weak constraints. In contrast, the tightest 
constraints come primordially, before collapse and heating occur, when the DM is highly 
non-relativisitic. In particular we will find that constraints on DM coupled to baryons at the 
time of recombination and DM coupled to baryons in protohalos can be very important, and 
this constraint will eliminate models whose charges are larger than about 10~ 6 , dependent 
on the mass of the DM. This constraint eliminates a broad class of models. 

On the other hand, this constraint does not eliminate DM with epsilon-charges that can 
give rise to a signal in direct detection experiments. We will find, however, that in the region 
where the DM could give rise to a signal in direct detection, one expects the DM to have been 
evacuated from the disk via supernova shock waves, and its re-entry to have been prevented 
by Galactic magnetic fields. Therefore, although direct detection experiments are extremely 



sensitive to small charges, we will find that charged DM, such as suggested in [20|, could not 
plausibly give rise to a signal in a direct detection experiment. We now go through these 
constraints in detail. 



III. RELIC DENSITY CONSTRAINTS 



We begin by discussing the constraints from the relic density. If the DM is non-thermally 
produced, its relic density depends on the production mechanism (for example, if the DM 
particle is produced via the decay of a mother particle, its relic density depends on the 
number density of the mother particle). In this scenario, constraints from the current relic 
abundance are highly model dependent. On the other hand, in the case of thermal relics, 
the DM density is simply determined by the thermally averaged annihilation cross section. 
As we will discuss explicitly, a charged DM consistent with all cosmological constraints must 
be non-thermally produced, unless it has additional interactions. We now review the relic 
density considerations. 

The DM can annihilate to photon pairs and to charged fermion pairs through the photon. 
In general, the charged DM can also carry other SM or hidden sector quantum numbers 
and annihilate through these channels as well. In our analysis we will not specify these 
additional interactions in detail; instead, we maintain a less model-dependent view. We 
assume DM is in thermal equilibrium in the early universe and require the irreducible 
annihilation processes not overly deplete DM. By considering only the annihilation channels 
induced by the electromagnetic charge of the DM, we can derive upper bounds on the charge 
e for a given mass mx- 

The annihilation cross sections of XX — > 77 and // at tree level are given by 



and 

ira 2 e 2 n / ml ( m 2 f , 



respectively, where qf is the charge of the SM fermion in units of electron charge and N c is 
the color multiplicity of the fermion. The total annihilation cross section of the DM particle 
at tree level is (cr an w re i) tot = (owz; re i) 77 + J2f (^W^rei)//- 

The tree level annihilation cross section is enhanced by the Sommerfeld effect in the low 



velocity dispersion 24j-|27j|. DM freeze out with Sommerfeld enhancement has been discussed 
in 28|, |29| . Since the mediator of the Sommerfeld enhancement is the standard model photon 
with zero mass, this enhancement never saturates. The enhancement factor for the tree level 
S-wave annihilation cross section is given by 

S= {aaa f n) l\, , (5) 

where v = v Te i/2 is the DM velocity in the center of mass frame. The thermally averaged 
total annihilation cross section including the Sommerfeld enhancement is given by 

Knfrcl) tot = (>an^rcl) tot 7^ J SV Tel e XXV "*' dv leh (6) 

where we assume a Maxwell-Boltzmann distribution for DM particle and xx = rnx/Tx with 
Tx as the DM temperature. The Sommerfeld enhanced annihilation tends to deplete DM 
particles with low velocity, which may distort the thermal distribution of the DM after kinetic 
decoupling. However, as we will show in the next section, the charged DM can couple to the 
thermal bath even during the recombination epoch, and therefore the Maxwell-Boltzmann 
distribution is a good approximation. 



Following the standard procedure to calculate the abundance of a thermal relic 30|, 131 
freeze out occurs when 

x f w m£--ln(mO (7) 



f = 0.038k(2 + K)m p imx{g/^)(a an v re i) 



tot J 



where x = mx/T, T is the temperature of the thermal bath, and g is the number of degrees 
of freedom of the DM particle; we take g — 4, for a Dirac particle. The value of k is chosen 
to match the numerical solution; we set K = 1. 

The present number density of the DM is the solution of the Boltzmann equation, which 
can be written as 

1 1 Fn~ f Xs {g*s/Jg7) (cTan^rel) 



Y(x s ) Y(x f ) V 45 



x 2 
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where Y = nx / s with s the entropy density, and x s = T s / nix with T s = 1 eV, where we stop 
the integration. Here we assume X and X have the same number density, nx = n x , and 
the total number density of the DM is their sum, nx+n x - Before kinetic decoupling, the 
DM temperature is the same as the thermal bath temperature and drops as a -1 , where a is 
the scale factor. After kinetic decoupling, the DM temperature drops as a" 2 . Therefore, the 
enhancement factor scales as S ~ x^ 1 and S ~ x~ 2 with respect to x before and after kinetic 
decoupling, respectively [28(. Since DM is tightly coupled to baryons through the massless 
photon until after the recombination epoch, S ~ x^ 1 over the entire range of integration of 
Eq. (Q. As we will show in Section IV, the elastic scattering rate rises as the temperature 
drops so that the DM may still be coupled to the thermal bath when the temperature is 
below 1 eV, even if e is chosen to satisfy the relic density bound. 



IV. STRUCTURE FORMATION AND CMB CONSTRAINTS 



A. Decoupling at the Recombination Epoch 



Charged DM particles interact with the Standard Model via a small coupling through 
the photon, so that Coulomb scattering can couple the charged DM to the baryon-photon 
plasma tightly even at low temperature. If the DM is still in kinetic equilibrium with the 
baryon-photon plasma during recombination, the DM density fluctuations can be washed 



out due to the radiation pressure and the photon diffusion (Silk damping 32j). The baryon 



acoustic peak structure will also be directly altered through the coupling. The effects of 



millicharged particles on CMB acoustic peaks have been discussed in Refs. |10L |33| . It was 
found that if the millicharged particles couple to the baryon-photon plasma tightly during the 
recombination epoch they behave like baryons, and CMB observations put an upper limit on 



their abundance fl m cph 2 < 0.007 (95%) [lOj. Here we assume that the DM is made of epsilon- 
charged particles and derive the relaxation time scale for the DM to reach kinetic equilibrium 
with baryons. To avoid damping effects on DM density fluctuations and CMB anisotropy 
constraints, we require that DM have completely decoupled from the photon-baryon plasma 
at the recombination epoch, and derive a bound on e for a given DM mass. 

We consider a DM particle that has momentum px in its comoving frame. After each 
scattering event the magnitude of the DM momentum changes by an amount 5px- The 
momentum transfer rate is thus given by 

d(^)/dt 

(Px) ' { } 
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Figure 1: Constraints from various sources, from top to bottom: (i) Scattering in the bullet 
cluster and NGC720, (ii) DM as a charged thermal relic, and (hi) DM virial processes, and (iv) 
recombination epoch. 

The thermally averaged momentum transfer per unit time is 

d(^)/dt =J2 Ub J d 3 VBdh x f(v B )f(v x )dn*^v rcl 6p 2 x , (11) 

b=e,p 

where daxb/dQ* is given by Eq. (El), rib is the number density of the baryon, and 5p x is the 
momentum transfer after one collision: 

^ = 2/^(1- cos 0*)- (12) 

Note that this quantity is reference frame independent. The thermally averaged momentum 
squared of the DM particle in its comoving frame is 

(Px) = J d 3 v x f(vx)(m x vx) 2 = ^m 2 x vl = 3m x T (13) 

for a DM particle in a thermal Maxwell distribution. To evaluate the thermal average for 
^reD we derive a general formula. For a given function of g(v Te {), we have 

/4 i _"?d + v v C A 

dv rel V* el ^— —je "06 ^oa+^ b > ob g(v Iel ), (14) 
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where we assume f(v a ,b) are Maxwellian distributions and t>o a ,b are the most probable 
velocities for the a and b particles, respectively. By using this general formula, we have 



d«)/dt = - n&QL t In (C m /2) , (15) 

b=e,p ( JL |- T 



mx m b 



where 8™ n is the cutoff in the integral. Its value is set by the maximum impact parameter 
due to Debye screening effects in the plasma. This maximum impact parameter is related to 
the minimum scattering angle through 



&ma x =^^COt(C72), (16) 
(Wei) 

where (jU&^rei) = 3T and cot(# mm */2) ~ 2/#™ m for small #™ in . The impact parameter of the 
scattering must not be larger than the Debye screening length, so we have b < 6 max = Xd- 
Thus 

or - (it) 



where A^ = a/T/ (4a cm irn e ) is the Debye length for the baryon plasma. Using Eq. (fl"3]) . we 
then have a momentum transfer rate of: 



i 

B=e, P 3m x Tl 



3TA 



D 



If the DM is tightly coupled to the baryon-photon plasma during recombination, DM density 
fluctuations will be damped. CMB observations also place strong constraints on the total 
abundance of charged particles in the tightly-coupled regime. Here we require that the 
relaxation time of the momentum transfer rate is larger than the Hubble time at the 
recombination epoch; i.e., 

T^Tr) > t R , (19) 

where tn ~ 3.8 x 10 5 years [34j, and Tr ~ 0.26 eV is the temperature at the recombination 
epoch. We take the baryon number density at recombination n e = n p = fLp c a R 3 /m p , where 
n b ~ 0.023/i" 2 , p c = 8.0992/i 2 x 1(T 47 GeV 4 , a R ~ 1/1091, and h ~ 0.71 Q- This constraint 
is plotted in Figured! In the above analysis, we have implicitly assumed that photons couple 
to baryons efficiently during the recombination epoch despite the presence of charged DM. 
We checked that electron-photon Compton scattering can keep baryons and DM particles in 
the kinetic equilibrium with photons. This is because the DM density is not far from the 
baryon density. The same Compton drag force also suppresses the growth of the DM density 
perturbations. 

Here we ignore the process of DM-photon Compton scattering. Since the cross section 
cr x 7 = 8 7ra em e4 /( m x) i s proportional to e 4 and the momentum transfer rate through this 
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process is also highly suppressed kinematically at low temperature, we expect that the bound 
derived from the DM-photon decoupling is weak. As shown in Ref. 35], CMB anisotropies 
and matter power spectrum requires ax-y/mx < 1CT 32 cm 2 GeV -1 . We can translate this 



limit to a bound on e as e < 0.49 (mj/(l GeV)) 3 / 4 , which is much weaker than the bound 
derived from the DM-baryon decoupling. 



B. Effect on the Dark Matter Virialization 



After recombination, radiation damping suppression is absent, but the efficient energy 
transfer between baryons and charged DM particles will modify the virialization process 
of the DM. Since baryons decouple from the thermal bath much later than DM particles, 
baryons are hotter than DM particles at redshift z ~ 30 when protohalos start to form. 
If there is a tight coupling between DM particles and baryons at this epoch, baryons will 
transfer energy to DM particles and heat them up. We can derive a bound on e by requiring 
the energy transfer time be longer than the DM virialization time. 

Eq. (jl8p is no longer valid for charged particles with slow motion in a neutral medium. At 
these late epochs, although it appears that the Born approximation condition e 2 a em /v Te \ < 1 
may still be satisfied due to the smallness of e, the charged DM particle typically has a 
wavelength larger than the Bohr radius of the hydrogen atom, and one must take into 
account the screening effect. This effect is analogous to the energy loss of a slow- moving 



ion in the neutral medium, a result first derived by Lindhard and Scharff [36J. Lindhard's 
approximation is valid when the impact parameter is larger than the Bohr radius. When the 
protohalo forms at redshifts z ~ 30, the DM velocity dispersion is (9(10~ 8 c). Its de Broglie 
wavelength is much larger than the Bohr radius under these conditions, and we expect that 
Lindhard's formula applies. We calculate the energy exchange of the charged DM in the 
hydrogen medium using Lindhard's formula, 

dEx Tin 



dl m e 



7r 2 e m x 
2.7183 mx + ran 



(20) 



where hh is the hydrogen number density, and we ignore the negligible effects of electron 



recoil 361 and other elements. The relaxation time scale is estimated as 



/ 1 dl \ 3 x 2.7l83m e (m H + m x ) T H T x , 01 , 
\v TC idE x / 4y/27i 3 ' 2 n H e V m H m x 

where Ex = mxv 2 el /2. We take \fT H / m H ~ 10~ 6 c, and \jT x /rnx ~ 10~ 8 c at z ~ 30. 

In the usual cold DM scenario, DM collapses and virializes at a redshift of z ~ 30. In 
over-dense regions the density is about 178 times larger than the average density at the same 



epoch 37j], and the violent relaxation time scale is 

rvir ~ (Gp tot )- 1/2 , (22) 
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where p to t = Px + Px an d Ptot ~ 178p to t = 178fixp c (l + z) 3 - Now we demand tx > t v[t and 
obtain an upper bound on e, 

£<2 , xl „-«(-±-), m 

which is shown in Fig. 1. 



V. DARK MATTER HALO CONSTRAINTS 



The strongest constraints on the coupling of DM to the photon come from scattering 
considerations rather than annihilations because of the large scattering cross-section 
enhancement at low velocities. In the previous section, we explored the effects at high 
redshift from observations of universe at recombination temperatures. The constraints are 
particularly strong in this regime. However, lower redshift observations can also be used to 
test the charged DM hypothesis. 



A. Elliptical galaxies 

If the DM self- interaction through Coulomb scattering is strong enough to create an (9(1) 
change in the momentum of DM particles within the age of galaxies, it will isotropize the 
velocity dispersion and lead to more spherical halos. The collisions also cause heat conduction 
from the hot outer parts to the cooler inner parts of DM halos, giving rise to the formation of 
a core with a shallow density profile. These expectations have been confirmed by simulations 



in the hard sphere scattering limit [38H43lj . In addition, observations of elliptical DM halos in 
clusters constrain self-interactions |44|. while observations of elliptical DM halos in galaxies 
provide the strongest constraints on self-interacting DM models , I45I . [46| . In this paper, 
we will follow the analysis of Ref. , I45I . [46| and use the ellipticity of NGC 720 to derive 
the upper bound of the electric charge of the DM. 

To estimate how the ellipticity of NGC 720 may be used to constrain the charge of the 
DM, we calculate the relaxation time due to momentum transfer. We then assume the 
relaxation time scale is the same as the time scale for isotropizing the mass distribution of 
the DM halo. By using Eq. (0Q), the thermally-averaged momentum transfer rate inside the 
halo can be evaluated from as 

r < = - 16 ™^ 4 3£r I^'^k'^^ r ' n/2) • (24) 

where p to t = Px + Px is the total DM density of the halo, and the minimal scattering angle 
Qmm j g gj ven by 0mm ^ 4. aeme / (3 mxV 2\ D } with X D = m x v q /(87re 2 a! em p t ot ) as the Debye 
screening length of the DM halo. 
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The elliptical galaxy NGC 720 is well-studied [47J, |48|. In Ref. [48[, X-ray isophotes were 
used to extract the ellipticity of the underlying matter distribution, and the DM halo of 
NGC 720 was found to be elliptical at 5 kpc and larger radii. At 5 kpc, the DM density is 



Ptot = 4 GeV/cm 3 , and the radial velocity dispersion t> 2 = t> 2 (r)/2 ~ (240 km/s) 2 45]. To 



derive the constraints on e for the given mx from the observed halo shapes, we require 

r; 1 > 10 10 years. (25) 

That is, the average time scale to create 0(1) change in the DM particle momentum must 
be greater than the galaxy's lifetime. This bound, weaker than the constraint derived from 
decoupling at the time of recombination, is depicted in Fig. [TJ 



B. The Bullet Cluster 



In the Bullet Cluster system, a subcluster has collided with and moved through a larger 
cluster. These clusters have three major components that each behave very differently during 
the collision: stars, gas, and DM. The visible stars pass through without colliding, but the 
highly collisional X-ray gas slows down significantly. Gravitational lensing shows that the 
DM tracks the stars, which are effectively collisionless. These observations have been used 
to place stringent bounds on the self- interaction of the DM [49]. These bounds are derived 
through different considerations including the offset between the gas and DM, the high 
velocity of the subcluster, and the survival of the subcluster after the collision. It turns 
out that the survival of the subcluster puts the strongest bound on the self-interaction of 

dm iif 



We follow the approach of Ref. (49J to derive bounds on the DM charge. The analysis of 
Ref. [49J is based on a hard sphere scattering cross section, but we relax this assumption. 
For Rutherford scattering, the subcluster experiences a net loss of DM particles if particles 
in both the main cluster and the subcluster have velocities larger than the escape velocity 
of the subcluster. We define the scattering angle 9 to be measured in the rest frame of the 
subcluster, which implies 9 = 6^/2, where 9* is the scattering angle in the frame of the center 
mass of two colliding particles. 

The particle loss condition detailed above can be satisfied if sin 9 is in the following range: 



^£ < s in#< Wl-% (26) 
Vi \j v( 

where vi ~ 4800 km/s is the velocity of the main cluster incoming particles before the 
collision and t> esc ~ 1200 km/s is the escape velocity of the subcluster. We assume that the 
subcluster sees the main cluster with a surface number density S m ~ 0.3 g/cm 3 0|, and 
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Figure 2: CoGeNT (blue), DAMA (green) allowed regions at 99% C.L. The CDMS-Si (yellow) line 
is included as a sample exclusion at 99% C.L. In the gray area, the charged DM is evacuated from 
the Galactic disk. Also shown the bound from Recombination epoch (red). Below the dotted line 
(black), charged DM may diffuse to the disk. 

demand that the particle loss fraction / be smaller than 30%, i.e. 

/ = £m fdO.^ = ^^44 (4- - ^ < 30%. (27) 
The Bullet Cluster bound is given in Fig. HJ 



VI. DIRECT DETECTION OF CHARGED DARK MATTER 



Because of the large enhancement of the scattering cross-section at low velocity, even DM 
with a very small charge can give rise to a large scattering cross-section in direct detection 



experiments. In [20J], it was found, for example, that a charge of e ~ 1CT 9 was sufficient 
to give rise to the relatively large signals in CoGeNT and DAMA. Thus, if correct, direct 
detection experiments have a potential to give rise to even tighter constraints on epsilon- 
charged DM with mass in the range mx ~ 10 GeV — 1 TeV. We find, however, that in the 
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range of charges where DM could give rise to a signal in a direct detection experiment the 
DM will necessarily have been efficiently evaporated from the Disk, and thus one expects no 
signal. We begin by a review of the signal in direct detection experiments. 



A. Direct Detection Basics 



The rate for scattering is 

dR „ p v f ,o „ , _ _ s da 



dh R m x J\$\ >Vmin dh R 



(2? 



where v min = vj^mEm ; anc l ^ N i s the reduced mass of the nucleus-DM system. We take the 
velocity distribution f(v,v e ) to be a modified Boltzmann distribution 

f(v, v e ) oc (e-^fht _ e -vlJv^ + (2Q) 

where explicit expressions for the velocity integrals from this distribution can be found 



m [50|. The additional term is to allow for a smooth cut-off of the velocity distribution 
near the Galactic escape velocity v esc . The Earth's speed relative to the Galactic halo is 
v e = v Q + Vorb cos 7 cos[o;(£ — to)} with v = t>o + 12 km/s, v or b = 30 km/s, cos 7 = 0.51, 
to = June 2nd and u = 27r/year. We take as a standard case Vq = 220 km/s, and we fix 
v esc = 500 km/s and the local DM density 0.3 GeV/cm 3 . 

A standard calculation relates the differential rate for scattering off nuclei to the scattering 
rate off a nucleus <ttv, 

da _ m N a N 

dE R ~ ' 1 J 

For the standard spin-independent case, this rate is related to a scattering off protons, a p , 
through 

iT _ u 2 N [f P Z + f n (A-Z)f 2 

Cat - Op — 7^ b V^r), (6L) 



where fi n is the DM-nucleon reduced mass and f p and are the DM couplings to the neutron 
and proton. We set f n = since the coupling is assumed to be through the photon. We make 
use of a Helm form factor F(E R ) = ^^e-^) 2 ^ 2 / 2 , where r = ((1.2A 1 / 3 ) 2 - 5s 2 ) 1/2 fm, 
with s — 1. The scattering cross-section off nuclei through the photon is 

167roi 2 e 2 Z 2 /4 

a N = ~ A , (32) 

which is to be inserted in Eq. ( 1281) to obtain the total rate as a function of energy. As an 
example of the typical DM charge e that can be probed with direct detection experiments, 
we show the constraints one obtains from the CDMS, DAMA and CoGeNT experiments in 
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Fig. (J2]). One can see that the viable region is well below the structure formation constraints 



labeled in the figure. Note that in Ref. 20], the mirror DM velocity dispersion depends 
on the particle mass, and is typically smaller than the rotation speed Vq = 220 km/s, so 
the allowed regions for both DAMA and CoGeNT as well as the excluded region for CDMS 
shift to larger DM mass (larger than 20 GeV) compared to the fitting presented in Fig. [2j 
Since the shift is more significant for the light target nuclei, the mirror DM model features 
a DAMA region which is not excluded by CDMS, as depicted in 20] . We next discuss how 
DM with charges in this range will have been evacuated from the disk at the present day, 
eliminating any possible signal in a direct detection experiment. 



B. Evacuation of Charged DM from the Disk 



The charged DM interacts with the magnetic fields of the Galaxy in addition to baryons 
in the disk. Since the large-scale magnetic field in the Milky Way is mostly parallel to the 
plane of the Galactic disk, the charged DM particle in the halo may not be able to penetrate 
the disk if its gyroradius is smaller than the height of the disk. This magnetic shielding effect 
for the millicharged particle has been discussed in the Ref. The gyroradius is given by 

where ~ 100 pc is the typical height of the Galactic disk. So the Galactic magnetic field 
prevents charged DM from entering the disk if e > 5.4 x 10 -13 (mx/GeV). 

It is possible that some quantity of charged DM can remain in the disk from the time when 
the disk formed. Subsequently, however, shock waves generated by supernova (SN) explosions 
can blow these particles out of the disk if the acceleration time scale (r acc ~ 10 7 years) is 
shorter than the cooling time scale (^]. The cooling time scale due to the scattering with 
electrons is given by 



m x m e v x 



In 



D 



(34) 



where we take n e ~ 0.025/cm 3 , the Coulomb logarithm In (/j, e v x \o/a em e) ~ 30 for the 
parameter range of interest. By demanding r coo i < r acc , we get e < 3.4 x 10~ 4 ^mx/GeV. 
Here we assume the epsilon- charged DM is efficiently accelerated by the Fermi mechanism. 
This is true when the gyroradius of the charged DM is smaller than the length of shock waves. 



Since the length of the shock wave can be ~ 100 pc 51J, as long as e > 5.4 x 10 13 (m x /GeV), 
the charged DM will be accelerated along with baryons over a time scale r acc ~ 10 7 years. 

Hence, if the DM charge is in the range 5.4 x 10~ 13 (m x /GeV) < e < 3.4 x 10~ 4 A/mx/GeV, 
the number density of the DM is negligible in the disk. We note that this constraint strongly 
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disfavors the charged DM explanations of DAMA and CoGeNT experiments, because the 
experimentally preferred value is e ~ 10~ 9 and nix ~ O(10 — 20) GeV. One possible way to 
relax this constraint is to consider the diffusion of DM into the disk. The Galactic magnetic 
field is not perfectly parallel, and in fact has a large nonperturbative turbulent component. 
The charged DM particles may diffuse into the disk as they interact with the turbulent 
magnetic field. The diffusion time scale can be estimated as t^s ~ 5Hl/(3R g vx) [52j. If 



the Tdiff is smaller than the acceleration time scale r acc , the charged DM may be able to 
diffuse to the vicinity of the earth and leave a signal in direct detection experiments. This 
signal is sensitive to the DM number density, which highly depends on the diffusion process 
and will in general be smaller than 0.3 GeV/cm 3 . We find that the charge has to satisfy 
e < 9 x 10~ 12 (nix /GeV) with vx = 270 km/s and B = 5 fiG, which is too small to fit 
DAMA and CoGent data. 



VII. CONCLUSIONS 



We have discussed the cosmological and direct detection constraints on charged DM. We 
considered in particular relic density, halo shape, large scale structure, recombination-era 
coupling, and direct detection constraints. We found that charged DM must have additional 
annihilation modes or be non-thermally produced if it is to satisfy the CMB constraints on 
DM couplings to baryons, which require the DM charge be smaller than e ~ 10 -6 for 1 
GeV DM (weakening to e ~ 10 -4 at 10 TeV). We discussed the possibility that one or more 
of these constraints is nullified by supernova shock waves blowing charged DM out of the 
disk, and showed as a result that DM with epsilon charge 10~ 9 cannot be an explanation for 
the CoGeNT or DAMA excesses, though the DM-baryon interaction cross-section is large 
enough. In addition, no signal in direct detection experiments could be expected in future 
experiments, as DM with large enough charge to generate a sizable DM-nucleus interaction 
cross-section would have been evacuated from the disk. 

While the idea of charged DM is in many ways an elegant one, its feasibility as a DM 
candidate, it appears, is strongly constrained. Further, these tight constraints also apply to 
any model where a massless dark photon kinetically mixes with the visible photon, so that 
fields charged under £/(1)em pick up a dark charge. This study presents constraints on a 
wide variety of hidden sector models that may be useful in the continued hunt for DM. 
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